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Abstract. We study the impact of modern stellar atmo- 
spheres that take into account the effects of stellar winds, 
departures from LTE and line blanketing ( " CoSta^ mod- 
els) on the ionization structure of H II regions. Results 
from a large grid of photoionization models are presented. 
Due to a flatter energy distribution in the He I continuum, 
compared to the widely used Kurucz models, generally 
higher ionic ratios are obtained. 

We find that N+/0+ and Ne++/0 ++ can be safely 
used as direct indicators of N/O and Ne/O abundance 
ratios in H II regions, over a wide range of astrophysical 
situations. 

The roughly constant observed value of Nc ++ /0 ++ 
ionic ratios in Galactic H n regions is naturally reproduced 
by photoionization models using CoStar fluxes, while Ku- 
rucz models at solar metallicity fail to reproduce this be- 
haviour. This gives support to ionizing fluxes from non- 
LTE atmospheres including stellar winds and line blan- 
keting. However, we also point out that tests of stellar 
atmosphere models from observations of H n regions are 
hampered by a lack of strong constraints on the ionization 
parameter. 

Key words: H n regions - ISM: abundances - stars: at- 
mospheres - stars: mass-loss - stars: early-type - ultravi- 
olet: stars 



1. Introduction 

Most of the recent work on H II regions modelling uses 
the Kurucz (1991) plane-parallel LTE model atmospheres. 
There are two reasons to this. One is that Kurucz's com- 
putations include line blanketing, which is indeed impor- 



tant in stellar atmospheres, especially at high metallici- 
ties. This effect was not included in earlier non-LTE stel- 
lar atmosphere modelling (e.g. Mihalas & Auer 1970). The 
other reason is that Kurucz models are available for a wide 
set of stellar temperatures and metallicities, which makes 
them extremely useful for any study aimed at reproducing 
the observed properties of H n regions in different environ- 
ments. 

However, it is recognized that Kurucz models cannot 
predict correctly the distribution of the radiation in the 
Lyman continuum of hot stars, not only because non- 
LTE effects are neglected, but also because the presence 
of winds in such stars is expected to have an impact on 
the outcoming stellar flux (Gabler et al. 1989, Najarro et 
al. 1996, Schaerer et al. 1996b). 

A new set of theoretical stellar spectra from the com- 
bined stellar structure and atmosphere models (referred 
to as CoStar models) of Schaerer et al. (1996a, b) has 
now been published (Schaerer & de Koter 1996, hereafter 
SdK96). The atmosphere, in particular, includes the ef- 
fects of stellar winds, departure from LTE as well as blan- 
keting. It is the purpose of this paper to explore the im- 
pact of these new ionizing fluxes on single star H II region 
models. 

In Section 2, we recall the main features of the CoStar 
models, and discuss the expected uncertainties. In Section 
3, we present a grid of photoionization models constructed 
for metallicities solar (Z©) and one tenth solar (Zq/10), 
using the CoStar models and the Kurucz ones, and dis- 
cuss the main differences. In Section 4, we present some 
implications on H II regions diagnostics. The possibility of 
testing these new stellar atmosphere models using obser- 
vations of H II regions is discussed in Sect. 5. Section 6 
summarizes the main conclusions of this study. 
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2. Comparison of CoStar and Kurucz model atmo- 
spheres 

2.1. Comparison of ionizing fluxes 

A detailed comparison of the ionizing fluxes from CoStar 
models with the plane parallel LTE models of Kurucz and 
also recent non-LTE models of Kunze et al. (1992) and 
Kunze (1994) is given in SdK96. Here we shall only briefly 
recall the most important differences between the CoStar 
and the Kurucz fluxes used in our study. Wind models 
similar to ours have also been presented by Sellmaier et 
al. (1996). 

As an illustration we show in Fig. 1 a comparison be- 
tween the CoStar and Kurucz model fluxes for two dwarf 
models with ~ 46000 and 35000 K respectively. Also 
indicated are the positions of some ionization potentials 
relevant to the discussion in forthcoming sections. For the 
CoStar model we plot the continous spectral energy dis- 
tribution with the spectral resolution used in the non-LTE 
calculation. Blanketing due to numerous spectral lines of 
all elements up to zinc is included in these models by 
means of an opacity sampling technique (using 20 A bins) 
as described in Schaerer & Schmutz (1994). The ioniza- 
tion edges of metals (mostly CNO, Ne, and Ar) found in 
other models (e.g. Kunze et al. 1992, Kunze 1994) are not 
treated in our models. On the other hand these may be 
overestimated if iron peak elements causing an important 
fraction of metal line blanketing arc not included. We note 
that the recent models of Sellmaier et al. (1996) also lack 
of pronounced metal ionization edges in the EUV, which 
seems to confirm our calculations. As discussed in SdK96, 
a more detailed analysis of this question will be possible 
with the inclusion of additional metals in the full non-LTE 
calculations (see e.g. de Koter et al. 1994, 1996). 

The largest differences between the Kurucz and CoStar 
models appear in the He n continuum above 54 eV: due to 
the velocity gradient in the stellar wind a strong non-LTE 
effect depopulates the He II groundstate and leads to a de- 
crease of the bound-free opacity above 54 eV (cf. Gabler 
et al. 1989). This explains the larger flux in the He II con- 
tinuum in non-LTE models accounting for stellar winds 
(CoStar) compared to the Kurucz models. Of more im- 
portance for the present study is the flux distribution at 
lower energies, where the ionization potentials of the most 
important elements observed in H II regions are located. 
The essential characteristics of the CoStar spectrum in the 
He I continuum (24.6 eV < E < 54 eV) are (cf. SdK96): 1) 
a natter spectral energy distribution, and 2) an increase 
of the total flux. The flattening of the spectrum mainly 
results from a wind effect whose strength decreases to- 
wards higher energies, since the continuum is formed at 
larger optical depths. In the considered temperature range 
non-LTE effects increase the total He I ionizing flux in 
both plane parallel models and those accounting for stel- 
lar winds. 



2.2. Adopted CoStar and Kurucz models 

The CoStar models are available for two sets of abun- 
dances: Solar (Z=0.020 in the notation of SdK96) and low 
metallicity (Z=0.004). We will compare these models with 
the Kurucz ones at solar and one tenth solar metallicity, 
which is the closest possible correspondance in metallicity. 

In the following, we will compare CoStar model atmo- 
spheres for main sequence stars (the ones labeled A2-F2 in 
SdK96) with the available Kurucz models having the clos- 
est effective temperature and gravity. We will also consider 
CoStar models for giant stars (the ones labeled A4-D4 in 
SdK96). 

The CoStar models provide not only the spectral dis- 
tribution of the ionizing radiation, but also the stellar 
radii. The radii of stars modelled with the Kurucz at- 
mospheres are obtained from the calibration of Vacca et 
al. (1996) for stars of luminosity class V and la, for main 
sequence and giants respectively. 

3. Photoionization models 

3.1. Definition of the grid 

The photoionization models are constructed with the code 
PHOTO, the latest description of which can be found in 
Stasihska & Leitherer (1996, hereafter SL96). A model 
computed with this version of PHOTO that can be di- 
rectly compared with the results obtained with other pho- 
toionization codes (presented at the Lexington workshop, 
Ferland et al. 1995) is also given in SL96 EJ. 

The ionization structure of a nebula is mainly deter- 
mined by the ionizing radiation field and by the nebu- 
lar geometry. In the simple case of a sphere of constant 
gas density n and constant filling factor e, the ionization 
structure of the heavy elements — for a given spectral dis- 
tribution of the ionizing radiation field — is essentially a 
function of the ionization parameter U = Qn/ (AirR 2 nc), 
where Qh is the number of stellar photons emitted per 
second above 13.6 eV, R is the external radius, and c is 
the speed of light. It can easily be shown that any combi- 
nation of Qh, n and e that keeps Qh ti e 2 constant results 
in a similar ionization structure. For each star, we have 
thus computed a series of 4 models of hydrogen density 
n = 10 cm~ 3 and e=l, in which the total stellar lumi- 
nosity has been multiplied by a factor 10 6 , 10 3 , 1, 10 -3 
respectively for the models numbered from 1 to 4. Doing 
this, one covers the entire range of expected ionization 

1 Extensive results (intensities of about 100 lines in the 
optical, UV and IR domains, mean ionic fractions and av- 
eraged electron temperatures) for the photoionization mod- 
els computed for this study are accessible by anonymous ftp 
on ftp.obspm.fr in the directory /pub/obs/grazyna/SS96 and 
most of them can be found as well on the CD-ROM accompa- 
nying the paper "A Database for Galaxy Evolution Modeling" , 
by Leitherer et al. (1996). 
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parameters, from the most compact H n regions (num- 
ber 1) to the most diluted interstellar medium (number 
4). Numbers 2-3 correspond roughly to the range of ion- 
ization parameters encountered in giant H II regions as 
found in SL96. The ionization structure of model number 
2 is the same as that of a model with a single ionizing star 
and a gas density of n = 10 4 cm~ 3 . From now on, models 
number 2 will be called the "reference models" □. 

Whithin the range of chemical compositions expected 
in H II regions, the ionization structure of a nebula is only 
little dependent of its actual chemical composition. 

In building our grid, the chemical composition of the 
gas was taken solar (0/H= 8.51 10~ 4 with the abun- 
dances of the remaining elements as in SL96) when dealing 
with solar metallicity stellar atmospheres, and one tenth 
solar(0/H=8.51 10 -5 ) when dealing with low metallicity 
atmospheres. 

3.2. Comparison of photoionization models using CoStar 
and Kurucz fluxes 

The comparison is made in a series of diagrams showing 
ratios of ionic fractions cc(X +n )/x(Y +m ) as a function of 
the stellar effective temperature B. Figures 2 and 3 rep- 
resent the solar abundance models using the CoStar fluxes 
and the Kurucz atmospheres respectively. Figures 4 and 5 
are equivalent to Figs. 2 and 3 respectively, but for the low 
metallicity models. In each figure, the "reference models" 
with different are connected with a thick continuous 
line for main sequence stars, and a thick dotted line for 
giants. For each star, models with different values of U are 
connected with a thin dotted line. 

In Figures 2-5 panels a - c give examples of ratios 
of ionic fractions that are supposed to indicate the stel- 
lar effective temperature (e.g. Rubin et al. 1995). Panel a 
shows x(He + )/x(H + ), which depends on the ratio of stel- 
lar photons above 24.6 eV and 13.6 eV. Panel b shows 
z(N ++ )/x(N + ), and Panel c shows ir(S ++ +)/x(S ++ ) 
for the solar abundance models and x(Ar +++ )/ir(Ar ++ ) 
for the low metallicity ones. While the x(He + )/x(H + ) 
is accessible from optical and radio measurements of 
H II regions (assuming a canonical helium abundance), 



Actually, a sphere may not be the best approximation to the 
geometry of H II regions. Therefore, we have also introduced 
a parameter /, which is the internal radius of the nebula, in 
units of the Stromgren radius for a full sphere at an electron 
temperature of 10000 K. We have computed series of models 
with / = 10~ 2 and / = 3. In the latter case, which is close to 
a plane parallel model, the average "excitation", as measured 
for example by [O m]/[0 n] tends to be lower than for a full 
sphere of the same U (i.e. with Qh, n, f and e combining into 
the same value of Qh n e 2 (1 + / 3 )~ 2 )- This effect is especially 
important for low effective temperature stars. In this paper, 
we present only results for / = 10~ 2 , but the interested reader 
will find the / = 3 models on the anonymous ftp account. 
3 a;(X +n ) is defined as: / n(X +n )n e dV / f n(X)/n c dV) 



cc(N++)/a;(N + ) and x(S +++ )/a;(S ++ ) are accessible from 
infra-red measurements. x(Ar +++ )/.x(Ar ++ ) can be ob- 
tained from optical lines, but only when the electron tem- 
perature is sufficiently high for the lines to be measur- 
able, that is to say only at low metallicities. Other ratios 
involving two subsequent ionization stages of the same 
element that can be compared in the same way with 
the observations are x(0 ++ )/x(0 + ), x(Nc ++ )/x(Nc + ), or 
x{S++)/x{S+). 

The fact that at a given T*, the CoStar models have 
higher fluxes in the ~ 30 eV - 54 eV region translates into 
higher ionic ratios in photoionization models computed 
with CoStar, as can be seen in Panels a - c of Figs. 2 
- 5. This is conspicuous in ratios involving ions of high 
ionization potential like S++ (34.83 eV) or Ar++ (40.74 
eV) . On the other hand, the effect is actually marginal on 
the x(He + )/x(H + ) ratio for main sequence stars, being 
noticeable only at effective temperatures smaller than ~ 
40000 K. The reason is that, at higher temperatures, there 
are enough photons above 24.6 eV to fully ionize helium, 
so that a;(He + )/a;(H + ) becomes independent of the model 
atmosphere. 

An effective temperature determined from an ionic ra- 
tio such as the ones mentioned above (assuming that one 
has some idea of the ionization parameter) would thus 
be lower using the CoStar models than using the Kurucz 
models. Note that a similar trend was found by Rubin 
et al. (1995) when they compared photoionization mod- 
els constructed using the non-LTE atmospheres of Kunze 
(1994) rather than the Kurucz ones. If indeed one were 
to use nebular models to determine the effective temper- 
ature, one would of course also have to take care that the 
volume sampled by the observations are comparable to the 
volume over which the mean ionic fractions are computed 
in the models. Furthermore one should also remember that 
the atomic data used in the photoionization model com- 
putations are still incomplete: dielectronic recombination 
is generally an important process in the ionization bal- 
ance, but the coefficients for ions of sulfur and argon are 
not available. Also, charge-transfer reaction rates are only 
approximately known for many ions. In spite of these dif- 
ficulties, it may be worthwhile noting that CoStar models 
predict a detectable [Ar iv] 4740 line in a wider range of 
effective temperatures than the Kurucz models. 

Another important issue is that with CoStar we have 
appropriate models for all the evolutionary stages of mas- 
sive stars. Models for giants tend to produce a higher ex- 
citation of the surrounding H n regions than what is ex- 
pected from extrapolating the dwarf models towards lower 
effective temperatures. This is clearly seen in Figs. 2 and 
4 where the curves representing the giants conspicuously 
depart from the curves corresponding to main sequence 
stars. The giant model with T + ~ 32000 K is a very strong 
case, where wind effects flatten the spectrum in the 20 - 
40 eV range (cf. SdK96). When using the available Kurucz 
models, the only difference between giants and main se- 
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quence stars is in the total luminosity, which induces only 
a negligible variation in the excitation of the surrounding 
nebulae, through an increase of the ionization parameter 
U. 

4. Some consequences for abundance determina- 
tions in H II regions 

Empirical abundance determinations in H II regions gen- 
erally use an ionization correction factor scheme first pro- 
posed by Peimbert & Costero (1969) based on the sim- 
ilarity of ionization potentials of several ions. For ex- 
ample, to derive the N/O ratio, one generally assumes 
that N/0=N + /0 + , and to derive the Ne/O ratio, one as- 
sumes that Ne/0=Ne ++ /0 ++ . Actually, photoionization 
models do not always agree with these empirical schemes 
(e.g. Rubin 1985, Stasiriska 1990). This is not surprising, 
since the ionization potentials of + and N + or Ne + are 
not exactly the same, and the various recombination pro- 
cesses do not have the same strength for these ions. Rubin 
et al. (1988) have strongly questioned the N/O ratios de- 
rived from optical measurements of H n regions (partly 
because of the implied ionization correction scheme) , and 
they consider the N/O ratios derived from far infrared line 
measurements assuming N/0=N ++ /0 ++ to be far more 
reliable. 

Panels d and e in Figs. 2-5 show x(N + )/x(0 + ) and 
x(N ++ )/a;(0 ++ ) respectively as a function of T*. We see 
that, actually, our models show the N + /0 + to be an ex- 
cellent approximation for N/O in all the range of effective 
temperatures considered, except for the Kurucz models 
at solar metallicity. In this case, the slope of the stellar 
spectrum aroud 30 eV is much larger than for the other 
models, inducing a differential effect between N + and + . 
Thus, if CoStar models describe the stellar radiation field 
more accurately than the Kurucz ones, the widely used 
empirical N/0=N + /0 + scheme remains a quite safe ap- 
proximation. Determining the N/O ratio using N ++ /0 ++ , 
on the other hand, gives rise to important biases since, at 
high effective temperatures, N ++ partly transforms into 
N +++ and at low effective temperature (when actually 
N + and + are the dominant ionic species), the predicted 
x(N ++ )/x(0 ++ ) becomes much larger than unity. 

Panel f in Figs. 2-5 plots x(Nc ++ )/x(0 ++ ) as a func- 
tion of T*. We see that, in all the "reference models", 
this ratio decreases with decreasing effective temperature, 
whether one uses the CoStar atmospheres or the Kurucz 
ones. However, the range in where x(Ne ++ ) / x(0 ++ ) is 
close to unity is much larger when using the CoStar at- 
mospheres, because they have flatter spectral energy dis- 
tributions in the 30 - 40 eV range than the Kurucz ones 
(see also next Section). Thus, if CoStar predicts correctly 
the stellar energy distribution in the Lyman continuum, 
one expects the Nc ++ /0 ++ to give directly the Ne/O ra- 
tios over a wide range of parameter space relevant to H II 
regions. 



5. Possible observational tests of stellar atmo- 
sphere models using H II regions 

Since the predictions on the ionization structure of H II 
regions depend so strongly on the type of model atmo- 
sphere, one might hope to use observations of H n regions 
to test the flux of these models in the Lyman continuum. 

Along this line, it has already been noted that the so- 
called [Ne ill] problem, which appeared when using the 
Kurucz model atmospheres, seemed resolved when using 
model atmospheres, which, similar to CoStar models, ac- 
count for stellar winds, and include a detailed non-LTE 
treatment and line blocking (Sellmaier et al. 1996, here- 
after SYPR96). This [Ne in] problem was made especially 
clear when Simpson et al. (1995; cf. references therein) ob- 
tained far infrared line measurements for a sample of H n 
regions and plotted the values of Ne ++ /0 ++ as a function 
of ++ /S ++ . They noted that the observations indicated 
an essentially flat Nc ++ /0 ++ in contradiction with the 
predictions of photoionization models using the Kurucz 
atmospheres (assuming no systematic variation in Ne/O 
or S/O). 

In Figures 6 - 9, we show the predictions of our model 
HII regions for Ne ++ /0 ++ to be compared directly to 
observations. Figures 6 and 7 represent the solar abun- 
dance models using the CoStar fluxes and the Kurucz 
atmospheres respectively. Figures 8 and 9 are equivalent 
to Figs. 6 and 7 respectively, but for the low metallic- 
ity models. The layout of the figures is similar to the 
one used in Figs. 2-5. Panel a plots Ne ++ /0 ++ = 
Ne/Ox(Ne++)/x(0++) as a function of 0++/S++ = 
0/Sx(0++)/x(S++). Here, Ne/O and O/S are the stan- 
dard abundance ratios used by Simpson et al. (1995) 
(Ne/O=0.2 and 0/S=47). In Figure 10, we plot the 
values of Ne ++ /0 ++ versus ++ /S ++ derived from in- 
frared line ratios for the HII regions observed by Simpson 
et al. (1995), using the atomic data referenced in SL96, 
and discarding measurements with upper limits only. Our 
models using CoStar have a similar behaviour to the ones 
of SYPR96), and are in much better agreement with the 
observations than models based on Kurucz atmospheres. 
Actually, for the hottest stars, we obtain slightly larger 
values of x(0 ++ )/x(S ++ ) at a given temperature T* than 
SYPR96. Indeed, the CoStar atmospheres seem to pro- 
duce somewhat higher fluxes in the 40 - 54 eV energy 
range than the atmospheres of SYPR96. Detailed com- 
parisons will be undertaken to understand the origin of 
these differences. It appears, however, difficult to distin- 
guish between both cases from nebular observations. 

Of course, one should bear in mind that different pho- 
toionization codes may produce slightly different ioniza- 
tion structures. A direct comparison can be made for the 
models using the Kurucz atmospheres. Our computations 
using PHOTO seem roughly in agreement with the ones 
of SYPR96, as can be judged by comparing our Fig. 7 
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with their Fig. 2, and are slightly different from the ones 
presented in Simpson et al. (1995). 

One can also perform a similar observational test of 
the model atmospheres by using optical lines, and plot- 
ting Nc ++ /0 ++ versus ++ /0 + , as shown in panel b of 
Figs. 6-9. This plot has the advantage of being indepen- 
dent of S/O. Figure 11 shows the observed values for a 
sample of Galactic H II regions (from Shaver et al. 1983, 
filled circles) and of Magellanic Cloud H II regions (from 
Pagel et al. 1978, open circles), derived using the atomic 
data listed in SL96. The trend seen in Fig. 11 is more com- 
patible with the CoStar models than with the Kurucz ones 
at solar metallicity. However, it must be noted that with 
0.25 solar metallicity models (more relevant to Magellanic 
Cloud nebulae), the problem with the Kurucz models is 
less severe than at solar metallicity. 

It is worthwhile recalling, that, in their abundance 
study of low metallicity H II galaxies from optical spec- 
troscopy, Masegosa et al. (1994) noted a completely dif- 
ferent [Ne in] problem. Their sample showed a slight ten- 
dency for the observed Ne ++ /0 ++ to increase as ++ /0 + 
decreases or as the H/3 equivalent width decreases. If 
Ne ++ /0 ++ were to measure the abundance ratio Ne/O, 
this would mean that Ne/O depends on the evolutionary 
stage of the H II galaxy. This is of course not acceptable, 
and Masegosa et al. suspected that the Ne ++ /0 ++ ra- 
tio was affected by charge exchange. Our models do take 
charge exchange into account. The fact that the charge 
transfer reaction ++ + H° is far more efficient than the 
Ne ++ + H° one produces a differential effect between oxy- 
gen and neon ions, which becomes more conspicuous as 
the proportion of neutral hydrogen becomes greater in the 
q++ region. This is seen in our models (Figs. 6-9): the 
Ne ++ /0 ++ ratio increases when the ionization parame- 
ter decreases. The effect found by Masegosa et al. (1994) 
is illustrated in Fig. 12, where we plot Ne ++ /0 ++ versus 
++ /0 + as derived in the sample of low metallicity H II 
galaxies studied by SL96 (which was based mainly on the 
Masegosa et al. line intensities). Because of these charge 
exchange reactions, it is thus difficult to test stellar model 
atmospheres from the diagrams shown above, if nothing 
is known about the ionization parameter. The only strong 
conclusion is that, for the Kurucz models to be consistent 
with the observations of nebulae with ++ /0 + < 1 or 
++ /S ++ < 10, one would require an unrealistically low 
ionization parameter, given the type of objects studied. 
On the other hand, atmosphere models including stellar 
winds, non-LTE effects and line blanketing naturally lead 
to approximately constant Ne ++ /0 ++ values over a large 
domain of stellar temperatures in agreement with obser- 
vations. 

Another difficulty in the testing of model atmospheres 
by using H n regions is that many observed H n regions 
are in fact complexes ionized by several O stars, up to 
hundreds of them. The best observational test would be 
in fact by modelling small single-star H n regions for which 



integrated line emission strenghts would be available. To 
our knowledge, such data do not exist presently, but will 
be available soon (L. Deharveng, private communication). 

6. Conclusions 

We have studied the impact of modern stellar atmospheres 
that take into account the effects of stellar winds, de- 
partures from LTE and line blanketing on the ionization 
structure of H n regions. Due to a flatter energy distribu- 
tion in the He I continuum, compared to the widely used 
Kurucz models, generally higher ionic ratios are obtained, 
especially for ions with high ionization potentials. For He 
the effect is noticeable for stellar effective temperatures 
below ~ 40000 K. 

Using the CoStar model predictions we find that 
N + /0 + and Ne ++ /0 ++ can be safely used as direct in- 
dicators of N/O and Ne/O abundance ratios, over a wide 
range of effective temperatures and ionization parameters. 

Comparisons of observed Ne ++ /0 ++ ionic ratios from 
IR line measurements and optical lines in Galactic H II re- 
gions (Shaver et al. 1983, Simpson et al. 1995) show that 
the CoStar models naturally reproduce the roughly con- 
stant observed value of Nc ++ /0 ++ , while Kurucz mod- 
els at solar metallicity fail to reproduce this behaviour. 
This gives a strong support to predictions of the ionizing 
fluxes from non-LTE atmosphere models including stellar 
winds and line blanketing, in agreement with the results 
obtained by Sellmaier et al. (1996). Compared to observa- 
tions of Magellanic Cloud H II regions (Pagel et al. 1978), 
the low metallicity Kurucz models, however, also show a 
reasonable agreement. 

We also point out that tests of stellar atmosphere mod- 
els from observations of H II regions may be inconclusive 
unless one has strong constraints on the ionization param- 
eter. Observations of integrated line strengths from single 
star H II regions with well known exciting sources should 
allow one to set much stronger constraints on the ionizing 
fluxes of O and B stars. This would be extremely helpful 
to asses the reliability of current models for O and B stars 
(cf. Cassinelli et al. 1995, Schaerer 1996). 
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Fig. 1. Comparison of the EUV continous spectral energy distribution from CoStar models (solid lines) and Kurucz models 
(dashed). Plotted is the photon number flux per unit energy as a function of the energy. Also shown are the ionization potentials 
of S+, S++, 0+, Nc+, and Ar++. Left panel: CoStar model D2: (T*,logg) = (46.1 kK, 4.05), Kurucz: (T*,logff) = (45 kK, 5.0). 
Right panel: CoStar model with (T*,logg) = (35 kK, 4.0) and remaining parameters from model Al, Kurucz: (T*,logg) = (35 
kK, 4.0) 
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Fig. 2. Predicted ionic fractions x(X +n )/x(Y +m ) as a function of the stellar effective temperature T* from photoionization 
models using CoStar atmospheres at solar metallicity (Z=0.020). At a given temperature the models 1 to 4, which differ by the 
ionization parameter, are connected by a thin dotted line. The "reference models" 2 are connected by the thick solid line for 
the main sequence stars and by the thick dotted line for giants respectively. See text for more explanations 





Fig. 3. Same as Fig. 2 for the Kurucz (1991) models at solar metallicity 
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Fig. 5. Same as Fig. 2 for the low metallicity (Z=0.002) Kurucz (1991) models 
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Fig. 6. Predicted ionic Ne ++ /0 ++ abundance ratio as a function of the ++ /S ++ ionic abundance ratio (panel a) and as a 
function of the ++ /0 + ionic fraction (panel b) for the solar metallicity CoStar models. Same symbols as in Fig. 2 




Fig. 7. Same as Fig. 6 for the Kurucz (1991) models at solar metallicity 
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Fig. 8. Same as Fig. 6 for the CoStar models at low metallicity (Z=0.004) 




Fig. 9. Same as Fig. 6 for the Kurucz (1991) models at low metallicity (Z=0.002) 
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Fig. 10. Observed ionic Ne ++ /0 ++ abundance ratios as a 
function of the ++ /S ++ ionic fractions in H II regions. Data 
from the IR line measurements of Simpson et al. (1995). To be 
compared with Figs. 6a - 9 a 
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Fig. 11. Observed ionic Ne ++ /0 ++ abundance ratios as a 
function of the 0++/0+ ionic fractions. Filled circles: Galac- 
tic H II regions, data from Shaver et al. (1983); open circles: 
Magellanic Cloud H II regions, data from Pagel et al. (1978). 
To be compared with Figs. 6b - 9 b 
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Fig. 12. Observed ionic Ne ++ /0 ++ abundance ratios as a 
function of the ++ /0 + ionic fractions. Data from the sample 
of low metallicity H II galaxies studied by SL96. To be com- 
pared with Figs. 6b - 9 b 
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